----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  F:/Dropbox (UBC-Umich)/Research_Projects/Project Javier and Natalia/TFP firm/TFP_firm_data_historical/dataverse_files/code/logs/figure_06_alt_transfer_costs.log
  log type:  text
 opened on:   7 Nov 2022, 11:50:58

. *-------------------------------------------------------------------------------
. 
. 
. use year isocode type D_A D_P using "${data}/estimates_agg_naics_sales_s1_base_withparent.dta", clear

. rename D_A D_A_new 

. rename D_P D_P_new 

. merge m:1 year isocode type using "${data}/estimates_agg_naics_sales_s4_base_withparent.dta", keepusing(D_A D_P)

    Result                           # of obs.
    -----------------------------------------
    not matched                             0
    matched                                27  (_merge==3)
    -----------------------------------------

. drop _merge 

. 
. *Assuming the technology transfer costs have an origin-specific component
. *-----------------------------------------------------------------
. clear all

. use year isocode type D_A using "${data}/estimates_agg_naics_sales_s1_base_woparent.dta", clear

. 
. 
. *Assuming the technology transfer costs have an detination-specific component
. *-----------------------------------------------------------------
. merge 1:1 year isocode type using "${data}/estimates_agg_naics_sales_s1_base_withparent.dta", keepusing(D_P)

    Result                           # of obs.
    -----------------------------------------
    not matched                             0
    matched                                27  (_merge==3)
    -----------------------------------------

. drop if _merge==2
(0 observations deleted)

. drop _merge

. 
. local grvar0 "gdp_pw_ppp" 

. local grvar1 "" 

. local rr ""

. local yy "2016"

. local yy_ref=`yy'       

. replace year=`yy_ref'
(0 real changes made)

. merge m:1 year isocode type using "${data}/aggregates_tfp_lp_klems_agg.dta" 
(note: variable year was float, now double to accommodate using data's values)
(note: variable isocode was str2, now str32 to accommodate using data's values)
(note: variable type was str5, now str9 to accommodate using data's values)

    Result                           # of obs.
    -----------------------------------------
    not matched                         2,985
        from master                         0  (_merge==1)
        from using                      2,985  (_merge==2)

    matched                                27  (_merge==3)
    -----------------------------------------

. drop if _merge==2
(2,985 observations deleted)

. drop _merge

. 
. foreach var in D_A  D_P {
  2. gen `var'_phi=(-1)*`var'*(0.2)
  3. *---------------------------------------
. }

. 
. local LHS1 D_A_phi

. local LHS2 D_P_phi

. local RHS ln_`grvar0'

.         
. sum `LHS1' `RHS'  `LHS2'

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
     D_A_phi |         27   -.2043066    .1622385  -.4886707   .0749743
ln_gdp_pw_~p |         27   -.3431099    .2740156  -.9310113          0
     D_P_phi |         27   -.1878562    .2014672  -.5504772   .1991635

. local ytitle "{&Delta}`=ustrunescape("\u03D5\u0303")'{sub:n}"

. local xtitle "{&Delta}y{sub:n}"

. local firm_emb_DA "{&Delta}D{sub:n}" 

. local firm_emb_DP "{&Delta}O{sub:n}"

. local title "Affiliates older than `edad' years"

. local title ""

. 
. reg `LHS1' `RHS'

      Source |       SS           df       MS      Number of obs   =        27
-------------+----------------------------------   F(1, 25)        =     12.07
       Model |  .222837263         1  .222837263   Prob > F        =    0.0019
    Residual |  .461517255        25   .01846069   R-squared       =    0.3256
-------------+----------------------------------   Adj R-squared   =    0.2986
       Total |  .684354518        26  .026321328   Root MSE        =    .13587

-------------------------------------------------------------------------------
      D_A_phi |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
ln_gdp_pw_ppp |   .3378563   .0972438     3.47   0.002     .1375789    .5381337
        _cons |  -.0883848   .0423907    -2.09   0.047    -.1756901   -.0010794
-------------------------------------------------------------------------------

. mat b = e(b)

. mat V = e(V)

. global c1: display %-03.2fc round(b[1,1],0.01)

. global se1: display %-03.2fc round(sqrt(V[1,1]),0.01)

. 
. reg `LHS2' `RHS' 

      Source |       SS           df       MS      Number of obs   =        27
-------------+----------------------------------   F(1, 25)        =     21.90
       Model |  .492780111         1  .492780111   Prob > F        =    0.0001
    Residual |  .562534708        25  .022501388   R-squared       =    0.4670
-------------+----------------------------------   Adj R-squared   =    0.4456
       Total |  1.05531482        26  .040589031   Root MSE        =       .15

-------------------------------------------------------------------------------
      D_P_phi |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
ln_gdp_pw_ppp |   .5024174   .1073601     4.68   0.000     .2813052    .7235296
        _cons |  -.0154718   .0468006    -0.33   0.744    -.1118595    .0809159
-------------------------------------------------------------------------------

. mat b = e(b)

. mat V = e(V)

. global c2: display %-03.2fc round(b[1,1],0.01)

. global se2: display %-03.2fc round(sqrt(V[1,1]),0.01)

. 
. display $c1
.34

. display $se1
.1

. display $c2
.5

. display $se2
.11

. 
. 
. global opt1 msymbol(circle_hollow) mlabel(isocode) mlabcolor(red)  msize(small) mlabsize(small)

. global opt2 msymbol(square_hollow) mlabel(isocode) mlabcolor(blue) msize(small) mlabsize(small)

. two (scatter `LHS1' `RHS', sort $opt1  mcolor(red)) (lfit `LHS1' `RHS', lcolor(red)) (line `RHS' `RHS', lcolor(none))  ///
>         (scatter `LHS2' `RHS', sort $opt2  mcolor(blue)) (lfit `LHS2' `RHS', lcolor(blue)) (line `RHS' `RHS', lcolor(none)),  ylabel(, format(%5.1f) labsize(medlarge)) xlabel(-
> 1(0.2)0, format(%5.1f) labsize(medlarge)) ///
>         title("`title'", size(small)) ytitle("`ytitle'", size(medlarge)) xtitle("`xtitle'", size(medlarge)) graphregion(fcolor(white) lcolor(white)) ///
>         legend(order(2 "`firm_emb_DA'=0: $c1 ($se1) (baseline)" 5 "`firm_emb_DP'=0: $c2 ($se2)") rows(2) size(medlarge) bplace(se) ring(0) region(lwidth(none))) 

. 
. graph export "${rmain}/fig_06.pdf", replace
(file F:/Dropbox (UBC-Umich)/Research_Projects/Project Javier and Natalia/TFP firm/TFP_firm_data_historical/dataverse_files/results/main/fig_06.pdf written in PDF format)

.         
.         
. log close 
      name:  <unnamed>
       log:  F:/Dropbox (UBC-Umich)/Research_Projects/Project Javier and Natalia/TFP firm/TFP_firm_data_historical/dataverse_files/code/logs/figure_06_alt_transfer_costs.log
  log type:  text
 closed on:   7 Nov 2022, 11:51:05
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
